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II i Abstract. Based on the non- linear logistic equation we study, in a qualitative 

Mh| and semi-quantitative way, the evolution with energy and saturation of the elastic 

differential cross-section in pp(j)p) collisions at high energy. Geometrical scaling occurs 
at the black disk limit, and scaling develops first for small values of the scaling variable 
|^|ftot.- Our prediction for da/ dt at LHC, with two zeros and a minimum at large |i| 
^ ' differs, as far as we know, from all existing ones. 
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Saturation phenomena are expected to dominate QCD physics at high energy and 
high matter density [H [2], as it may happen at LHC and cosmic rays at ultra high 
energies. Non hnear differential equations include, in a natural way, saturation effects. 
This happens with the well known logistic equation[3], which can be seen as a simplified 
version of the B-K equation[l]. See[5] and [6] for discussions on evolution and saturation. 

We shall concentrate here in the evolution of the imaginary part of the impact 
parameter elastic amplitude, or the profile function r(6, R) = Im B[b, R), where b is the 
impact parameter, related to angular momentum i by 

b^l=i, (1) 

where ^/s is the centre of mass energy, and R is an increasing with energy radial scale 
parameter. Partial wave unitarity constrains r{b,R): 

0<T{b,R)<l. (2) 

We now write two logistic equations, in R and b, respectively: 

and 

where 7 > is pratically a constant. From (jS]) one sees that dT/dR > and that, for 
fixed b and F > 0, F reaches the black disk limit: F = F^ = 1. FromQ one sees that 
in general F is a decreasing function of b and that, for large b, F decreases, as expected, 
exponentially (F ~ exp —b/'j), saturation occurring first at small b. 
A solution of ([3]) and (j4]), not the most general one, is: 

mR) = — ^ir-T- (5) 

exp ^^— ^ + 1 
The total and elastic cross-section are written as 

atods) = 2t, jv{b,R)db'' (6) 

and, neglecting real part contributions, 

crez.(s) = vr f \V{b,R)\''db\ (7) 



respectively. The imaginary part of the elastic amplitude Im F(t, R) is the Fourier-Bessel 
transform of F(6, R) and the differential elastic cross-section is written as 

da_ ^ al^ \mF{t,Rf 

dt 167rImF(0,/?)2' '^ '' 

It should be noticed that in regions of energy where 7/-R — const., as it happens 
at ISR energies (20 < y/s < 60 GeV), T{b,R) satisfies geometrical scaling [7], 

nb,R) ,77=^ r(/3), (9) 
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and, 



nb,R) j^np)^{: -._ , (12) 



P^^, (10) 

^^R^f{tR^))\ (11) 

As 7, contrary to R, does not show in general a consistent dependence on energy, 
in the hmit i? — )■ oo, •y/R — )■ and one obtains again scahng [8], 

1, 0</3< 1 

^0, /3 > 1 

with atot. ~ CTei. ~ -R^ and aei./atot. -^ const. = 1/2. 

It should be also noticed that the parameter i? in ([5]) separates the region of negative 
curvature, b < R, from the region of positive curvature, b > R. In fact, R plays the role 
of the angular momentum L, used in the proof of Froissart bound [.9| by Martin and 
collaborators [10], that separates the region that contributes in a significant way to the 
total cross-section, i < L, from the region that is negligible, C. > L: 

^=4^ = ^lnf-), (13) 

with y/t^ = 2mT^. 

When comparing our model to experimental da/dt one finds that 7 takes values 
always of the order 7 ~ 1 mb ' , while i? is a monotonically increasing function of 
energy. 

In fact 7, the parameter that controls the low density region, b > R, can be seen 
as a measure of the range of the interaction, with 7 ~ (2m^)~^ in the Yukawa picture. 
The evolution of the cross-sections with the energy is controlled by the single parameter 
R{s), the effective impact parameter radius. 

When studying the evolution of the amplitude ([5]) with energy one finds three 
regimes (see Fig. [1]): 

i) v^<20GeV, Fig. [Ha). 

This is the region corresponding to linear evolution, with F small and with 
exponential behavior, da/dt being a monotonically decreasing function of —t. 

ii) 20 < ^ < 63 GeV (ISR energies). Figs. [Hb) and c). 

In this region a dip, which is a minimum, not a zero, appears at —t ^ 1.4 GeV^ 
and slowly moves to the left as energy increases. Conventional wisdom says that the 
dip results from a zero: interference between one-Pomeron and two-Pomeron exchanges 

m- 

iii) 500 < v^ < 1.8 TeV, Figs. Ed), e) and f). 

In this region the minimum becomes negative, originating a pair of zeros. Instead 
of the clean second maximum of region ii) one has now a kind of shoulder, but with a 
cross-section higher by an order of magnitude. 
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Figure 1. da/dt as function of — < at different energies, showing the sequence: no 
structure in a), one minimum in b) and c), two zeros in d) and e) and two zeros and 
one minimum in f) at LHC. Values of 7 and J?: a) 7 = 1.020, R = 1.972; b) 7 = 1.026, 
R = 2.171; c) 7 = 1.090, R = 2.259; d) 7 == 1.000, R = 2.575; e) 7 = 1.016, R = 2.496; 
f) 7 = 1.078, R = 2.683; g) 7 = 1, i?. = 3.770. Data from [H]. Dashed line: only 
imaginary part contribution. Full line: the real part of the amplitude is included. 



In Fig. [Hg) we have also included our expectation for LHC (assuming atot. ~ 110 mb 
see [I3] for expected range of values - and 7 = 1 mb ' ). Our LHC curve clearly shows 
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Figure 2. d{a/dt) / {da/dt{{))) as a function of tlie scaling variable |i|CTtot. showing 
the approach to black disk geometrical scaling from small to larger values of the scaling 
variable. Scaling only applies to the imaginary part of the amplitude. The parameter 
7 was put equal 1. 



how the evolution towards the black disk continues: from a monotonically decreasing 
curve at large —t a minimum starts developing which at some stage generates a pair 
of zeros to join the previous pair. And so on! In the black disk limit we just have a 
sequence of pairs of zeros. 

At high energy, y^ > 60 GeV, when app ~ app is not difficult to use the derivative 
dispersion relations p^ to estimate the real part contribution to the differential cross 
section. In Figs. [T]c) to [T]g) we show, in full line, da/dt including the real part 
correction. The real part contribution is important at the zeros of F. 

In Fig. |2]we show the geometrical scaling plot (see ( ITTl) ) of ^/^{t = 0) as function 
of the variable |t|o"iot. for different values of R and for 7 = 1 mb^'^, and the black disk 
limit. The way the approach to the scaling curve is achieved seems clear: as the energy, 
or R, increases scaling is satisfied for larger |t|crtot. values. The LHC curve corresponds 
to atot. = 110 mb, aei./atot. - 0.28. 

In Fig. |3] we show the correlation between (Tei.l<^tot. and atot. in comparison with 
data. Note the transient geometrical scaling at ISR energies, aei./atot. — const.. 

It should be mentioned that the curves shown in Fig. [1] are not fits to data, but 
represent quahtative descriptions. This means that a^i. and atot.i which are controlled 
by the first points in the low | t \ region, are not constrained. That is the reason why 
the values of R{s) in Fig. [T]do not coincide with the values used in Fig. [31 For instance, 
at ISR energies, in Fig. [1], cxez. and atot. are larger by a factor of the order of 15% relative 
to the true values (in Fig. [3|). 

Recently, in several papers [15l [161 EZ] the questions of soft physics were addressed 
from different points of view. However, contrary to the present paper, no attempts 
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Figure 3. aei./crtot. as a function of atot. malcing use of ([S]). Data points, from left to 
right: ^i = 23.5 GeV, 7 = 1.065, i? = 1.710; ^/s = 30.6 GeV, 7 = 1.090, R = 1.731 
Vi = 44.9 GeV, 7 = 1.095, i? == 1.805; ^/s = 52.8 GeV, 7 = 1.107, R = 1.810 
Vi = 200 GeV, 7 = 1.123, R = 2.100; ^/i = 540 GeV, 7 = 1.080, i? = 2.498 
Vs = 900 GeV, 7 = 1.026, i? = 2.662; ^ = 1800 GeV, 7 = 1.046, R = 3.050. We 
expect for LHC (Tei./<Jtot. = 0.28 and atot. = HO mb. Data from [11]. 



were made to explain the qualitative aspects of da/dt evolution with energy. Instead, 
interesting problems related to inelastic diffraction, as Higgs production, were addressed. 
The black disk perspective varied significantly from [15] to [IE] : the black disk limit being 
already present at LHC energies in [15], but occuring at extremely high energies in [T6] . 
For a discussion on the black disk limit see also [T8] . 

The relevance or not of the Froissart bound (see [19] and [20]) and a dynamical 
interpretation of it (see [21] and [22]) are still matters open to discussion. 

Finally, we summarize our work. Starting from the non-linear logistic equation we 
obtained a solution for the high-energy imaginary part of the amplitude, and we were 
able to describe in a qualitative and semi-quantitative way the essential features of the 
evolution of the differential elastic cross-section with energy, namely the sequence: no 
structure in |t|, one minimum, two zeros and so on. Our prediction for da/dt at LHC 
energies is different from all the ones we are aware of (see, for instance, [23]). 
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